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Abstract 

We evaluate the effective interactions in a ffuid of electrons moving in a plane, on the 
approach to the quantum phase transition from the paramagnetic to the fully spin-polarized 
phase that has been reported from Quantum Monte Carlo runs. We use the approach 
of Kukkonen and Overhauser to treat exchange and correlations under close constraints 
imposed by sum rules. We show that, as the paramagnetic fluid approaches the phase 
transition, the effective interactions at low momenta develop an attractive region between 
parallel-spin electrons and a corresponding repulsive region for antiparallel-spin electron 
pairs. A connection with the Hubbard model is made and used to estimate the magnetic 
energy gap and hence the temperature at which the phase transition may become observable 
with varying electron density in a semiconductor quantum well. 
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Electronic fluids with an essentially two-dimensional (2D) dynamics, as are realized in 
semiconductor heterostructures, present a very rich phenomenology especially at low density, 
where the interactions between carriers come to play a crucial role [1]. This fact continues 
to motivate theoretical studies of the simplified system consisting of a gas of electrons de- 
scribed by the basic Coulomb Hamiltonian with a uniform neutralizing background (EG). 
More generally, this has been for many decades a basic reference model for calculations of 
electronic structure in molecular and condensed- matter systems [2]. 

Quantum Monte Carlo (QMC) simulation studies have been highlighting the formation 
of spontaneous spin polarization in the EG at strong coupling. A continuous transition from 
the paramagnetic to spin-polarized states has been reported to take place in the 3D EG with 
increasing coupling strength [3, 4], before a first-order transition into a ferromagnetic Wigner 
electron crystal occurs. Spontaneous partial spin polarization ("weak ferromagnetism" ) has 
been observed in doped hexaborides [5]. In the 2D EG, on the other hand, the quantum 
simulation studies indicate a first-order transition to a fully spin-polarized ("ferromagnetic") 
fluid before crystallization, no evidence being found for partially spin-polarized states [6 - 
8]. 

Theoretical studies of the equation of state described by the ground-state energy E{n, Q 
of the EG with spin densities and as a function of the electron density n = n-f + n ^ and 
of the spin polarization ( = (nj — n±)/n are clearly called for in both dimensionalities. The 
ground-state energy is in turn determined by the pair distribution functions, stressing the 
importance of two-body electron-electron scattering via effective interactions in the many- 
body system [9, 10]. 

Exchange between parallel-spin electrons and correlations from the Coulomb repulsions 
are crucial in determining the pair distribution functions. These effects are embodied in the 
so-called local-field factors (LFF) entering the expressions of the charge and spin susceptibil- 
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ities of the EG [11]. The LFF play the same role as vertex corrections in accounting for the 
difference between the effective potentials experienced by an electron and their mean-field 
value. From their knowledge one can also calculate the quasi-particle self-energies and the 
effective electron- electron interactions, as is well known from work on the unpolarized EG 



In this letter we are concerned with the effective electron-electron interactions in the 
paramagnetic phase of the 2D EG as functions of the coupling strength on the approach to 
the magnetic quantum phase transition. We adopt the linear-response approach of Kukko- 
nen and Overhauser [9] to construct these interactions from recent determinations of the 
LFF [13, 14], embodying their asymptotic behaviours from exact sum rules and information 
from QMC studies. The difference between the jj and jj effective interactions is related 
to the magnetic susceptibility and in the low scattering-momentum region most strikingly 
reflects the approaching phase transition in its behaviour with increasing coupling strength. 
After pointing out how these results could be used to reveal the magnetic phase transition 
in a theoretical calculation based on the solution of Schrodinger equations for two-electron 
problems, we establish an approximate connection with a simple mean-field solution of the 
Hubbard model for itinerant electrons and use it to estimate the temperature of the mag- 
netic transition in the 2D EG at the density predicted by the QMC studies. 

Following, therefore, Kukkonen and Overhauser [9], the effective interaction potential 
between two electrons with spin indices <j\ and a 2 is written in momentum space as 



where by convention the product <Ji • a 2 equals +1 for parallel spins and -1 for antiparallel 
spins. The functions Vo(q) and J(q) are expressed in terms of the bare Coulomb potential 
v(q) = 27ie 2 /q and of the ideal-gas Lindhard susceptibility Xo(<?) as 



[12]. 



v aia 2 (q) = V (q) - J{q)a 1 ■ a 2 



(1) 



Vo(q) = v{q) 



l + [l-G + (q)]G + (q)v(q)xo(q) 
1 + [1 - G + (q)}v(q) X o(q) 



(2) 
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and 

J{q)=v{q)- — 7T (3) 

1 - G-(q)v(q)x {q) 

Here, G+(q) and G-(q) are the charge-charge and spin-spin LFF in the paramagnetic 
2D EG. It is important to remark for later use that the denominator in Eq. (3) also enters 
the microscopic magnetic susceptibility x(q) °f the EG, which is given by 

Uq) = ^ (4) 

X[q) l-G-(q)v(q)xo(q) K) 

in units such that the Bohr magneton equals unity. 



The limiting behaviours of the LFF at low and high momenta are determined by a set 
of exact sum rules [15, 16]. Considering first G+(q), these are the compressibility sum rule 
yielding 



HmG +(? ) 



7T 8V2 



d 2 e c {r s 
dr 2 



1 de c (r s 



kp 



(5) 



with the Fermi momentum given by kp = \j2ixn in terms of the areal density n , and the 
Kimball-Niklasson-Holas relation yielding 

r s d[r s e c (r s )} q 



lim G+(q) = [1 - g(0)\ r- , , . 



(6) 



In these equations e c (r s ) is the correlation energy per particle as a function of the density 
parameter r s = (ima^) -1 / 2 with a# the effective Bohr radius, and g(0) is the pair distribution 
function at contact. The corresponding limiting values of G-(q) are determined by the 
magnetic susceptibility sum rule, 



limGL(g) 



d 2 e r (r x 



C=o. 



Q_ 

kp 



7i 2V2 dC 
and by the Santoro-Giuliani-Holas relation, 

lim G4 q )=9(0)-^ d[rs£ / rs)] ^-. 



(7) 



(8) 



Ample data from QMC studies are available for the correlation energy of the 2D EG as a 
function of density [7, 17] and for its long- wavelength magnetic susceptibility determining 
the coefficient in Eq. (7) [8, 18 - 20]. The value of g(0) is also known over a wide range 
of density [21], from the solution of the two-body scattering problem and from many-body 
calculations in the ladder approximation. 

Complete analytical expressions for both LFF have been obtained in the density range 
r s < 10 [13, 14], using all the above information and further QMC data from [18 - 20] to 
describe their behaviour in the region of intermediate wave number. We have now extended 
this work to determine the LFF at r s = 20, with the results that are shown in Figure 1. Es- 
pecially relevant for this purpose have been the newly available QMC data on the magnetic 
susceptibility [8], which reaches the value x(0)/Xo(0) — 25 at r s = 20. The corresponding 
values of the effective interactions between parallel-spin and antiparallel-spin electrons at 
several values of r s are shown in Figures 2 and 3, respectively. Some smoothing of apparently 
noisy features has been made in obtaining these results. 

As was shown in Ref. [13, 14] and as again assumed in constructing G-(q) in Figure 1, 
the long- wavelength behaviour reported for this LFF in Eq. (7) holds over a rather wide 
range of scattering momentum q, extending almost up to q = 2kp. The plateau shown by 
the effective interactions in this range of momentum is, therefore, essentially determined by 
the value of the magnetic susceptibility x(0)- As already noted, this thermodynamic param- 
eter is approaching very large values as the paramagnetic phase of the 2D EG approaches 
the quantum phase transition to a fully spin-polarized state, that the latest QMC evidence 
places at r s m 25 [8]. The effective interaction between parallel-spin electrons in Figure 

2 becomes correspondingly strongly attractive in the momentum range q < 2k F , while a 
repulsive structure grows at momenta above 2kp. In parallel to this behaviour, and as is 
evident from Eq. (1), the effective interaction between antiparallel-spin electrons in Figure 

3 becomes strongly repulsive in the momentum range q < 2kp. Although the quantitative 
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details of these results are rather sensitive to the precise value of x(0), their significance 
is evidently correct: the approach to the transition into an ordered phase having a space- 
independent order parameter is driven by the emergence of effective interactions favouring 
spin alignments over a scale of distances which is becoming long compared with {2kp)~ l ■ 

In Figure 2 we have also reported the values of V^(q) as calculated in the fully spin- 
polarized state at r s = 40 [22], to show that the effective interactions regain a "standard 
shape" after the magnetic phase transition has taken place. That is, the electron-electron 
interactions in momentum space are repulsive and show a rapid variation as q goes through 
2k F , corresponding to Friedel oscillations in coordinate space. 

As was shown in the work of Takada [23] (see also Richardson and Ashcroft [24]), the 
effective electron-electron interactions may be used to assess the possibility of superconduc- 
tive s-wave or p-wave pairing in the EG through the solution of an Eliashberg equation in 
the absence of phonons. It is natural to ask whether they may play a similarly useful role 
in regard to a magnetic phase transition. In fact, Ortiz et al. have stressed in their work [4] 
that a precise signal of the transition of an electron fluid to a magnetically ordered state is 
carried by the spin- spin radial distribution function, which is proportional to the difference 
<7Tf(r) — <7fj.(r) between the radial distribution functions for the two spin populations. Near 
a magnetic transition the function g^(r) — gn(r) can be expected to change sign at short 
range, going from negative in the deep paramagnetic phase (where it reflects a preference for 
spin alternation) to positive as magnetically ordered domains emerge. Methods to evaluate 
pair distribution functions in a many-body system through the solution of a Schrodinger 
equation for a two-body scattering problems have drawn some attention in the recent lit- 
erature [10, 21, 25] and could easily be adapted to the study of magnetic phase transitions 
with the help of the effective potentials that we have determined in this work. 

Here we adopt a much simpler and cruder route in order to estimate the magnetic 
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energy gap A which would be required for spin flips at zero momentum transfer in the 
magnetically ordered phase. This estimate is based on the remark that Eq. (4) for the 
magnetic susceptibility takes at low momenta the form that is obtained in the Hubbard 
model for a system of itinerant electrons [26]. The value of the Hubbard parameter U 
appropriate for the 2D EG near the magnetic transition can therefore be estimated as 



that is, taking r s w 25 and x(0)/xo(0) ~ 25 as from Ref. [8], A « 3.KT 3 Ryd and 
T c ~ 400 K for the EG model. This value is reduced by about three orders of magni- 
tude when one takes account of an effective mass correction and of screening by a back- 
ground dielectric constant in a semiconducting material. For instance, for the high-mobility 
GaAs/GaAlAs heterostructures studied by Yoon et al.[27] we estimate T c rs 1 K. 

In summary, we have evaluated the effective interactions in momentum space between 
electrons in the paramagnetic phase of the 2D EG at various densities and shown how their 
shape reflects the emergence of magnetically ordered domains on the approach to the quan- 
tum phase transition into a fully spin-polarized phase. We have also given an estimate for 
the transition temperature at the density predicted by the QMC runs. 

This work was partially funded by MIUR under the PRIN2001 Initiative. 



U « Jim [G-(k)v(k)\ = Xo \0) ~ X-\0) 



(9) 



Setting A rs nU for a fully spin-polarized state we estimate 




(10) 
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Figure caption 




FIG. 1. Charge-charge and spin-spin local field factors, G+(q) and G-(q), as functions of g/fep 
in the 2D EG in the paramagnetic state at r s = 20. The straight dashed lines show the asymptotic 
behaviours determined by sum rules and exact relations (Eqs. (5) - (8)), while the full line show 
an approximate interpolation suggested by earlier studies at lower r s [13, 14]. 
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FIG. 2. Effective interaction potential in momentum space for parallel-spin electrons in the 2D 
EG in the paramagnetic state at various values of r s from 5 to 20. The dots report the effective 
interaction in the fully spin-polarized 2D EG at r s = 40 [22]. 




FIG. 3. Effective interaction potential in momentum space for antiparallel-spin electrons in the 
2D EG in the paramagnetic state at various values of r s from 5 to 20. 
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